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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  were  obtained  by  ARO, 
Inc.,  AEDC  Division  (a  Sverdrup  Corporation  Company),  operating  contractor  for  the 
AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project  Numbers 
E32C-04  and  E32C-04A.  The  Air  Force  project  manager  was  Mr.  Eules  L.  Hively.  The 
manuscript  was  submitted  for  publication  on  September  28,  1979. 
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1.0  INTRODUCTION 

In  order  that  measurements  ean  be  made  of  structural  and  aeromechanical  turbine 
engine  behavior,  it  is  imperative  that  there  exist  a  turbine  engine  test  facility  capable  of 
simulating  flight-maneuver-mduced  gyroscopic  and  inertial  loadings  on  operating  gas 
turbine  engines  (Refs.  1  through  4).  The  requirements  for  such  testing  have  been  defined, 
and  a  test  facility  has  been  described  that  would  meet  these  objectives.  A  preliminary 
design  concept  of  this  facility,  called  the  Turbine  Engine  Loads  Simulator  Facility 
(TELS).  is  shown  in  Fig.  I.  The  facility  is  a  large,  rotating  arm  with  a  variable-position 
engine  mounting  structure  and  counterweight.  The  arm  is  of  a  truss-like  structural 
configuration  and  is  about  40  ft  long  from  the  engine  centerline  (most  extended 
position)  to  the  center  of  rotation.  It  is  designed  to  rotate  at  angular  speeds  of  up  to  3.5 
radians/sec.  The  engine  support  structure  also  supports  and  positions  the  X-ray  radiograph 
monitoring  system,  which  evaluates  engine  structural  deflections  or  motions  during 
imposed-loadings. 

To  operate  properly,  test  engines  must  ingest  air  of  adequate  flow  quality.  Air 
entering  the  engine  compressor  and  fan  must  have  acceptably  low  total-pressure  and 
total-temperature  spatial  distortions  or  nonuniformities.  For  these  conditions  to  exist,  the 
engine  inlet  must  operate  unstalled  or  unseparated,  and  the  jet  thrust  deflector  must 
direct  hot  exhaust  gases  away  from  the  engine  inlet.  These  conditions  must  be 
incorporated  into  the  design  of  the  TELS  facility  so  that  routine  loads  testing  can  be 
satisfactorily  performed  with  the  engine  running  reasonably  dose  to  actual  flight 
maneuvers  operating  points. 

The  objective  of  the  present  investigation  and  analysis  was  to  assess  inlet  flow 
quality  for  turbine  engine  operation  in  TELS  by  evaluating  inlet  and  jet  deflector 
performance  in  crossflows.  Another  objective  was  to  devise  a  method  for  predicting 
engine  performance  subject  to  the  predicted  inlet  flow  quality. 

2.0  ANALYSIS 


2.1  FLOW  SIMULATION 

A  survey  of  crossflow/cross  wind  effects  on  gas  turbine  engine  inlets  established  that 
an  extensive  and  thorough  study  of  inlet  behavior  is  being  performed  by  researchers  at 
the  NASA  Lewis  Research  Center,  Cleveland,  Ohio.  The  main  thrust  of  the  NASA 
research  lies  in  the  design  and  development  of  inlets  tor  V/STOL  aircraft  applications. 
One  important  product  of  this  research  program  has  been  the  development  of  technology 
lor  predicting  three-dimensional  potential  flow  in  plane  or  in  axisymmetric  inlets  subject 
to  uniform  but  arbitrarily  directed  crossflow.  This  technology  and  its  applications  for 
interpreting  inlet  behavior  in  cross  winds/crossflows  are  described  in  Refs.  5  through  11. 
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Another  aspect  of  this  technology  for  inlet  evaluation  is  the  development  and  use  of  a 
modern  finite-difference  boundary-layer  solution  procedures  (Ref.  12)  to  identify 
flow-separation  events  in  inlets.  Thus,  once  the  potential  main  or  bulk  inlet  flow  field  has 
been  found,  then  the  viscous  behavior  of  the  flow  in  inlets  can  be  predicted,  based  on 
boundary-layer  theory. 

The  technology  necessary  to  evaluate  inlet  behavior  in  crossflows  was  obtained  from 
N.  O.  Stockman  of  NASA  Lewis  Research  Center.  The  system  was  made  operational  at 
AEDC  by  G.  W.  Lewis,  ARO,  Inc.  In  addition,  Lewis  added  computer-plot  subroutines 
for  displaying  longitudinal  and  traverse  contours  of  constant  velocity  and  static  pressure. 
A  set  of  inlet  designs  was  selected  for  evaluation;  these  designs  covered  the  normal  range 
of  inlet  configurations  currently  available  for  crosswind  engine  testing.  Three  geometries 
from  this  set  are  shown  in  Fig.  2.  The  compressor  face  was  a  model  of  the  FI 00  engine 
configuration.  The  flow  fields  were  calculated  for  the  worst  case  of  crossflow,  =140 
ft/sec.  This  cross-flow  maximum  velocity  corresponds  to  the  tip  velocity  of  the  TELS 
arm  at  an  angular  velocity  of  3.5  radians/sec.  The  angle  of  attack  of  the  inlet  for  the 
worst-case  condition  corresponds  to  90  deg,  the  vertical  engine  position  shown  in  Fig.  1. 

The  calculations  provided  (1)  the  spatial  location  of  the  stagnation  lines  in  and 
around  the  inlets,  (2)  the  general,  three-dimensional  velocity  field  over  the  inlet 
configuration,  and  (3)  the  flow  inside  the  inlet  up  to  the  simulated  compressor  face 
station.  A  boundary-layer  analysis  of  the  potential  flow  field  revealed  unsatisfactory 
performance  of  the  first  two  inlet  designs,  A  and  B  (Fig.  2).  The  flow  separated  from  the 
forward  or  upstream  lip  of  these  geometries  in  the  meridional  or  symmetry  plane.  For 
the  third  geometry  design,  C  (Fig.  2),  the  boundary-layer  analysis  indicated  no  flow 
separation  or  stall.  Hence,  this  wide,  flared-bellmouth  design  would  operate  unstalled  for 
the  worst  case  of  TELS  crossflows.  The  free-stream,  inlet  streamline  pattern  is  shown 
schematically  in  Fig.  3.  Figure  3  shows  a  section  of  the  three-dimensional  stream  tubes,  a 
view  in  the  plane  of  the  cross-flow  velocity.  The  indicated  stagnation  points  are  points  on 
stagnation  streamlines.  In  the  meridional  plane  the  rear  stagnation  point  lies  on  the  inner 
edge  of  the  flared  lip.  This  may  indicate  a  possible  problem  with  flow  stability; 
calculations  made  with  the  computer  simulation  showed  this  stagnation  point  location  to 
be  sensitive  to  variations  in  engine  airflow  rates  and  crossflow  velocity. 

This  flared-bellmouth  geometry  is  only  one  example  of  the  geometries  that  would 
produce  acceptable  flow  quality,  and  it  may  not  necessarily  represent  the  best  design.  An 
optimization  of  the  inlet  design  may  be  performed  in  future  studies  using  these 
computational  tools.  Measures  of  quality  to  be  optimized  are  discussed  in  Section  3.0. 
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The  qualitative  nature  of  the  unstalled  velocity  field  in  the  flared  inlet  is  shown  in 
Fig.  4.  The  calculated  potential  flow  field  was  interpolated  to  provide  contour  plots  of 
constant  total  velocity  magnitude,  axial  velocity  component,  and  radial  velocity 
component.  The  velocity  field  is  not  symmetrical  about  the  centerline,  indicating  that  a 
distorted  static-pressure  field  will  exit  in  the  inlet.  However,  data  exist  showing  that  this 
may  not  adversely  affect  engine  performance. 

In  Ref.  13,  a  turbofan  engine  with  an  asymmetric  inlet  was  tested  at  high  crossflow 
velocities  (up  to  160  knots)  and  targe  angles  of  attack  (up  to  120  deg).  The  engine 
operated  satisfactorily  without  flow  separation  at  conditions  even  severer  than  the 
worst  ease  of  TELS.  Though  asymmetric,  the  inlet  looked  quite  similar  to  inlet  B  in  Fig. 
2,  but  it  bad  a  thicker  tip  facing  the  crossflow  direction.  Thus,  it  appears  that  the 
presence  of  the  engine  can  considerably  reduce  inlet  separation,  compared  to  analytic 
predictions,  and  can  handle  distorted  velocity  in  static  pressure  fields  reasonably  well. 

2.2  GAS  INGESTION 

A  mathematical  simulation  of  the  turbojet  engine  and  jet  deflector  was  formulated 
to  estimate  the  severity  of  the  recirculating  exhaust  gas  ingestion  problem.  The  simulation 
was  developed  as  a  two-dimensional,  finite-difference  model  of  a  turbojet  engine  and  an 
exhaust  gas  deflector.  The  computational  mesh  encompassed  the  surrounding  flow'  field 
far  enough  upstream  and  downstream  of  the  engine  to  simulate  the  undisturbed  far-field 
upstream  and  a  “fully  developed”  downstream  condition.  (However,  an  actual 
three-dimensional  flow  will  not  reach  a  fully  developed  state  but  dispersed  in  some 
manner  depending  on  the  shifting  of  the  wind  over  the  existing  terrain.)  The  fully 
developed  flow  condition  is  a  useful  mathematical  boundary  condition.  The  undisturbed 
free  stream  (far  above  the  engine)  and  the  ground  plane  formed  the  other  boundaries  of 
the  simulation.  Figure  5  gives  the  model,  the  mesh  spacing  system,  and  the  lengths 
incorporated  in  the  simulation. 

The  crossflow  interaction  with  the  deflected  exhaust  jet  and  the  engine  inlet  flow 
was  computed  in  this  mesh  system  by  a  finite-difference  solution  of  the  full,  elliptic, 
Navier-Stokes  and  energy  equations  subject  to  imposed  boundary  conditions.  In  this 
study,  the  effects  of  compressibility  were  neglected.  However,  the  solution  procedure 
written  included  a  Poisson-type  equation  for  the  static  pressure  that  may  be  solved 
simultaneously  with  the  other  field  equations.  The  finite-difference  solution  procedure 
was  used  in  a  Limited  parametric  study  of  various  effects  on  the  severity  of  engine  inlet 
ingestion  of  hot  exhaust  gases.  The  present  study  was  limited  to  an  investigation  of  the 
effects  of  cross  wand  velocity  at  a  constant  turbulent  viscosity  level  on  gas  ingestion  for  a 
fixed-geometry,  flat-plate  jet  deflector.  Although  the  deflector  geometry  was  flat,  curved 
deflectors  can  be  simulated  using  different  boundary  conditions  along  the  deflector 
surface. 
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The  velocity  and  temperature  distributions  shown  in  Fig.  6  are  typical  results  of  the 
present  numerical  investigation.  In  this  case  a  sine  function  was  used  to  describe  the 
deflector  boundary  conditions.  (The  actual  distribution  is  provided  in  Sections  A6.0  and 
A7.0  of  the  appendix).  These  results  show  no  significant  gas  ingestion  problem  in  the 
engine  inlet  plane.  The  temperature  of  the  gas  entering  the  inlet  is  only  slightly  higher 
(10°R)  than  the  free-stream  temperature  of  and  at  a  uniform  spatial  profile.  Figure  7 
shows  isotherms  in  the  flow  in  the  engine  tailpipe  vicinity.  The  flat  plate  is  assumed  to 
rise  to  the  exhaust  gas  temperature;  thus,  the  flow  above  the  deflector  is  hot.  Obviously, 
sensitive  instrumentation  should  not  be  placed  aft  of  or  in  the  wake  of  the  deflector. 

The  value  of  the  Reynolds  number  selected  for  this  flow  was  much  lower  than  that 
normally  used  for  turbulent  mixing  based  on  simple  jet  mixing  theory.  The  lower  number 
was  used  in  order  to  account  for  the  turbulence  created  by  the  rotation  of  the  TELS 
arm.  The  arm  can  be  expected  to  cause  a  local,  well-mixed  crossflow  of  extremely  high 
turbulence  and  hence  high  effective  viscosity.  Calculations  were  made  with  values  of 
turbulent  Reynolds  number  an  order-of-magnitude  higher  than  those  displayed  in  Figs.  6 
and  7,  and,  though  the  velocity  and  the  temperature  fields  changed,  the  inlet  temperature 
profiles  remained  essentially  flat.  Similar  results  were  obtained  for  calculations  performed 
for  crossflow  velocities  ranging  from  25  to  140  ft/sec. 

The  analysis  carried  out  in  this  study  suggests  that  with  a  careful  design  of  engine 
mounting  structures,  there  is  little  chance  of  exhaust  gas  ingestion  for  test  cases  that  use 
flat-plate-type  deflectors  located  behind  a  vertically  mounted  engine  in  a  uniform 
crossflow  velocity.  The  analysis  provides  a  framework  for  investigation  of  other 
parameters  affecting  gas  ingestion  problems.  Future  work  will  require  investigation  of  the 
following  problem  areas: 

1.  engine  orientation  in  crossflow, 

2.  deflector  geometry  (i.e.,  such  geometries  as  bucket,  angled  plate,  or 
splitter  plate), 

3.  turbulence  transport  model  (i.e.,  in  which  either  one-  or  twoequation 
models  can  be  incorporated  in  the  solution  procedure), 

4.  effect  of  exhaust  gas  deflector  proximity  on  the  exhaust  jet  static  pressure 
field  in  engine  nozzle  exit, 

5.  buoyancy  effects  caused  by  gas  temperature  variations, 

6.  three-dimensionality  effects,  including  the  scavenging  or  pumping  effect  of 
the  rotating  TELS  arm,  and 

7.  engine  exhaust  conditions. 
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The  present  study  neglected  the  bouyancy  effects  in  the  Navier-Stokes  equations. 
Tliis  allowed  the  analysis  to  provide  a  “worst-case”  description  of  the  flow  field  since,  in 
actual  TELS  operation,  bouyancy  effects  would  act  to  lift  the  hot  gases  away  from  the 
engine  inlet.  The  appropriate  source  terms  are  indicated  in  the  Navier-Stokes  equations 
(Appendix  A);  these  terms  can  be  easily  incorporated  into  the  solution  procedure. 

As  remarked  before,  there  is  a  simple  method  for  simulating  the  effects  of  other 
deflector  geometries  by  distributing  the  velocity  profiles  along  the  flat-plate  surface.  For 
example,  sine  or  cosine  distributions  can  be  defined,  with  the  u  and  v  components  equal 
to  zero  at  the  deflector  plate  centerline,  but  increasing  to  maximum  negative  or  positive 
values  at  the  ends  of  the  deflector.  This  procedure  artificially  modifies  the  resultant 
exhaust  jet/deflector  flow  field.  A  velocity  prediction  representing  a  given  deflection 
geometry  may  be  obtained  by  first  predicting  -  with  a  method  like  that  of  Ref.  14  — 
the  velocity  profiles  of  the  essentially  inviscid  flow  into  and  out  of  a  given  deflector  and 
by  then  using  the  exit  profiles  as  the  mathematical  boundary'  condition  in  a  program  like 
the  one  developed  for  the  present  study.  Figure  8  illustrates  this  method. 

2.3  DISTORTION  INDICES  AND  ENGINE  MATHEMATICAL  MODELS 

Work  has  been  done  for  several  years  on  the  formulation  and  interpretation  of  gas 

turbine  engine  inlet  airflow  quality.  Parameters  have  been  developed  whose  values 

correlate  with  measurable  changes  in  engine  behavior.  These  indices  are  assigned  values  by 
use  of  mathematical  procedures  to  evaluate  engine  inlet  flow  quality.  The  flow  quality  is 
defined  theoretically  by: 

1.  the  circumferential  and  radial  patterns  of  steady-state  total  pressure, 

2.  the  circumferential  and  radial  patterns  of  steady-state  total  temperature, 

> 

3.  swirl  strength  and  distribution. 

4.  average  one-dimensional  total  pressure  flustrations  of  “high”  (greater  than 
one-per-compressor  rotor  revolution)  and  “low”  frequency,  and 

5.  average  one-dimensional  total  temperature  fluctuations. 

The  flow  quality  has  not  been  related  either  to  the  velocity  ot  to  inlet  static 
pressure  fields,  but  experimentation  and  a  variety  of  rationale  have  produced  several 
distortion  index  parameters.  A  future  standard  method  for  the  industry  is  described  in 
the  Society  of  Automotive  Engineers’  Aerospace  Recommended  Practice  (ARP)  1420, 
Ref.  15.  The  mathematical  procedures  outlined  in  this  standard  are  recommended  for 
future  calibration  or  calculation  of  engine  distortion  indices.  However,  each  current 
engine  may  have  indices  calculated  with  procedures  unique  to  that  engine.  For  example, 
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Table  1  lists  the  engine  mathematical  models  available  at  AEDC  as  recently  surveyed  by 
the  author.  Of  these  mathematical  models,  only  those  for  the  FI 00  and  FI 01  engines 
could  account  for  distortion  effects.  The  FI  00  engine  was  selected  as  representative  of 
the  engine  type  tested  in  TELS. 

The  purpose  of  the  present  study  was  not  to  review  these  indices  but  to  learn 
whether  such  indices  existed  in  mathematical  models  of  engines  likely  to  be  tested  in 
TELS  and  to  determine  a  method  for  predicting  or  interpreting  engine  performance  in 
TELS  tests.  In  carrying  out  this  research  it  was  found  that,  in  addition  to  the  correlation 
methods  for  gaging  engine  performance,  both  experimental  and  theoretical  studies  are  in 
progress  to  formulate  methods  for  predicting  engine  response  to  pressure  and  temperature 
perturbations.  Compressor  models  are  being  formulated  and  evaluated  to  gain  insight  into 
compressor  response  to  inlet  distortion  (Refs.  16  and  17). 

By  using  such  models  to  predict  engine  performance,  it  may  be  possible  to 
accurately  gage  the  engine  response  to  inlet  flow  distortions  before  TELS  testing. 
Therefore,  basically  two  methods  were  found  for  estimating  the  F100  engine  performance 
in  TELS.  One  method  is  based  on  the  parallel  compressor  model  for  the  F100  described 
by  Walter  and  Shaw  (Ref.  17).  The  other  method  requires  using  two  computer  programs. 
First,  the  Pratt  and  Whitney  deck  CCD  1087  is  used  to  convert  input  total  pressure  and 
temperature  profiles  to  a  set  of  inlet  distortion  indices  or  K-factors.  Then,  these  K-factors 
are  input  to  the  Pratt  and  Whitney  FI 00-3  engine  performance  program,  deck 
CCD-l  116-3.0,  to  define  the  resultant  performance.  Similar  methods  are  available  or  can 
be  developed  for  other  current  turbine  engines. 

Some  cautionary  remarks  are  needed  to  qualify  the  foregoing  conclusions.  First,  data 
exist  that  indicate  that  inlet  performance  can  be  improved  with  an  attached  engine.  The 
engine  has  an  upstream  effect  on  the  flow,  not  only  delaying  or  minimizing  flow 
separation  but  also  modifying  the  velocity  and  static  pressure  field  in  the  inlet.  Therfore, 
analytical  predictions  of  inlet  performance  must  be  cautiously  applied  when  interpreting 
flows  or  defining  distortion  indices  for  inlets  attached  to  operating  turbojet  engines. 
Second,  the  effects  of  temperature  nonuniformity  are  twofold.  Temperature 
nonuniformity  may  affect  the  aerodynamic  behavior  of  the  flow  through  the  compressor 
by  coupling  with  the  pressure  field  to  produce  aerodynamic  stall,  or,  in  variable-cycle 
engines,  a  localized  temperature  may  affect  sensors  that  establish  blade  or  stator  angles  of 
attack,  thereby  causing  engine  component  mismatch.  Third,  an  area  of  concern  for  TELS 
is  the  interaction  of  the  engine  with  the  inlet  flow  during  engine  power-up.  Procedures 
must  be  devised  for  accelerating  the  engine  while  bringing  the  TELS  facility  to  speed, 
without  creating  transient  engine-damaging  How  conditions  (i.e.,  those  conditions  that 
exist  during  transient  operating  procedures  that  could  damage  the  engine,  as  opposed  to 
the  undamaging,  final  steady-state  conditions). 
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3.0  CONCLUSIONS 

An  investigation  was  performed  to  assess  inlet  flow  quality  for  turbine  engine 
operation  in  THLS.  First,  it  was  confirmed  that  an  axisymmetric,  flared  bellmoutli 
operates  satisfactorily,  i.e.,  unseparated.  Second,  a  theoretical  study  of  the  exhaust  gas 
ingestion  in  TELS  indicates  that,  for  the  simple  engine/blasl  deflector  arrangement 
simulated  in  the  study,  the  gas  ingestion  problem  is  insignificant.  Third,  a  method  has 
been  devised  lor  evaluating  the  performance  of  engines  to  be  tested  in  TELS.  The 

method  is  based  on  the  calibrated  mathematical  engine  models  that  incorporate  distortion 
indices. 
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Figure  1.  Preliminary  engineering  concept  of  TELS 
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Figure  3.  Approximate  streamline  pattern  into  flared  bellmouth 
inlet,  no  inlet  separation. 
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a.  Lines  of  total  velocity  magnitude 
Figure  4.  Lines  of  isovelocity  in  meridional  cross  section  of  flared 
inlet  showing  nonuniformity  of  velocity  field. 


b.  Lines  of  axial  velocity  component 
Figure  4.  Continued. 


c.  Lines  of  radial  velocity  component 
Figure  4.  Concluded. 
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Turbulent  Reynolds  Number 
Rt  -  2.2  Based  on  b  =  0.5  ft 

(see  Appendix) 

Crossflow  Velocity:  25  ft/sec 


a.  Profiles  of  axial  velocity  upstream  of  engine 
Figure  6.  Theoretical  flow  fields  around  the  TELS 
engine-exhaust  deflector  system. 
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b.  Axial  distribution  of  vertical  velocity  along 
two  horizontal  planes 

Figure  6.  Continued. 
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Figure  7.  Isotherms  around  the  engine-exhaust  deflector  system. 
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Figure  8.  Models  of  exhaust  gas  deflectors. 
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Table  1.  Engine  Mathematical  Models  at  AEDC 


Engine  Designation 

In  Current  Use 

J57-P-21/23 

Yes 

J75-P-1 7 

Yes 

YJ101-GE-100 

No 

TF30-P-3 

Yes 

TF30-P-7 

No 

TF30-P-100 

Yes 

TF33-P-7A 

Yes 

TF41-A-1 

Yes 

TF34-GE-100 

No 

F101-GE-100* 

Yes 

J85-GE-21 

No 

FI 00-PW-1 00* 

Yes 

J79-GE-15 

Yes 

TF39-GE-X 

No 

FI 07-WR-1 00 

Yes 

*Indicates  mathematical  models  known  to 
include  distortion  indices. 
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APPENDIX  A 

COMPUTATIONS  OF  EXHAUST  GAS 
INGESTION  BY  ENGINES  IN  TELS 

A1.0  DESCRIPTION  OF  NUMERICAL  SIMULATION  OF  GEOMETRY 


The  geometry  simulated  is  a  rectangular  engine  with  a  flat-plate  jet  deflector  in  a 
rectangular  flow  region.  The  engine  is  3  ft  wide  and  1 5  ft  long.  The  6-ft-wide  jet 
deflector  is  positioned  5  ft  behind  the  engine.  The  flow  region  is  80  ft  long  and  55  ft 
high.  The  boundary  conditions  necessary  for  the  computation  are 


Jet  inlet  velocity 
Jet  exhaust  velocity 
Jet  exhaust  temperature 
Crossflow  velocity 


400  fps 
1,900  fps 
1 ,500°  R 
25  to  140  fps 


The  computational  mesh  is  uniform  and  square  with  0.5-ft  spacings. 

A2.0  SYSTEM  OF  GOVERNING  EQUATIONS 


These  equations  are  for  turbulent,  incompressible,  heated  flow's.  The  equations  are 
the  Navier-Stokes  equations  and  the  energy  equation  in  plane.  Cartesian  coordinates. 


A2.1  NAVIER-STOKES  EQUATIONS 


du  du  I  dp 

u  —  +  v—  = - — i 

r}\  ay  p  n\ 


(A2.1) 


dv  dv 

U  —  +  V  — 
dx  dy 


p  dy 


“g  0  'P  ~PJ 


A2.2  ENERGY  EQUATION 


dr  <?t 

U  -  +  V  - 

dx  dy 


kE  d2T  d2T 
pCp  dx2  "*  dy2 


where  the  dissipation  function  is 


<t>T  = 


(A2.2) 


(A2.3) 


(A2.4) 


29 


A  EDC-TR  -79-83 


A2.3  PRESSURE  EQUATION 


For  the  pressure  field,  the  Poisson’s  equation 


V2p  = 


(A2.5) 


given  by  Roache  (Ref.  18)  is  solved.  The  density  is  related  to  the  local  pressure  and 
temperature  by 


P  =  P  RT 

A3.0  STANDARD  FORM  OF  THE  GOVERNING  EQUATIONS 

The  equations  are  all  rearranged  in  the  form 


,  d<t>  ,  .  _  n 

2^  —  blrf,~  b2(j>  dv  ~  c4  “  0 


^x2  24>  i?v2  ^  ^ 


The  coefficients  are  listed  below  for  each  variable' 


(A2.6) 


aK 

$ 

ro 

\ 

V 

$ 

s. 

u 

i 

i 

u/«E 

v/vE 

-(1 /p) 9p/9x 

V 

i 

i 

U/Vj. 

v/vE 

-(1/p)9p/3y 

T 

i 

i 

U?r  /VE 
rE  E 

'V'’* 

VPV*T 

P 

i 

i 

0 

0 

<t> 

These  equations  are  solved  by  the  Chien  “decay  function”  method  of  finite 
differences  (Ref.  I  9).  The  continuity  equation  is  not  explicitly  included  in  the  governing 
set  of  equations;  however,  in  the  solution  procedure,  global  checks  on  mass  conservation 
are  performed. 


A4.0  FINITE-DIFFERENCE  SOLUTION  PROCEDURE 


The  finite-difference  method  is  point-by-point,  Gauss-Siedel  Relaxation  of  a  system 
of  algebraic  equations  of  the  form 


<h 

‘•j 


l.J 


C. 

i  - 


.  * 

l.j  i- 


C  ,  O 
i,j+  l  i 


,]+  l 


(A2.7) 


j-i 


l-j 
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where  C,  ,  are  the  coefficients  for  the  following  five-point  computational  molecule 
i  +  l,j  =  Eastern  Node 


C  ,  .  =  C,  =  2 

At,.  A\  -  It |  U\t 

m  t-  1  -J  G 

CR 

i  —  1 , j  =  Western  Node 

C  i  —  ] -j  =  CW  = 

(A]  /  AX  -  1)1 .  BX) 

CR 

i,j  +  1  =  Northern  Node 

C  .  =  Cv.  = 

i.j -  1  X 

(A2.  AY  -  13-2  BY) 

cn 

i,j  —1  =  Southern  Node 

c.,-,  =  cs  = 

(A2/AY  H  B2/BY) 

CR 

where  the  A’s  and  B’s  are 

A1  =  ai<t, 

A2  =  a, 

B1  =  % 

B2  =  b0 

AX  =  GX  DX2 

AY  =  GY  DY 

BX  =  2  DX 

BY  =  2  DY 

Finally,  the  CR  parameter  is  calculated  from  the  following: 

CR  -  (2.0/ A XI  -  <2.0/ AY) 

The  point-wise  source  term,  Dij,  is  as  follows: 

- \Jn 

The  decay  functions,  GX  and  GY,  are  defined  according  to  Chien  (Ref.  19): 

(  1  -0.0625  Rx2  if  |  Rx  |  <  2 

GX  =  \ _ l _ L_  if  |  Rx  |  >  2 

(  |  Rx |  Rx2 


(ATS) 

(A2.9) 

(A2.10) 

(A2.ll) 


(A2.12) 

(A2.13) 

(A2.14) 
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(  1  -  0.0625  Hy2  if 

GY  -  ]  2  ,  ,f 

C  |  Hy  I  Ry2 


Rx  and  Ry  are  cell  Reynolds  numbers  defined  as 


Rx 


a 


ft 


Kyi  <2 
Ry  |>2 


(A2.I5) 


(A2.16) 


%  A>' 


2$ 


A5.0  TURBULENCE  TRANSPORT  THEORY 

v\ i  is  an  effective  turbulent  viscosity  given  by  Schlichting  (Ref.  20). 

"E  =  k//VJETRJET  =  0.0285  VJET  Rjet  (A2.17) 


where  VjEt  and  Rjet  are  the  turbojet  exhaust  velocity  and  jet  radius. 


PrF  is  the  effective  turbulent  Prandtl  number  defined  as 


P 

rE 


(A2.18) 


For  most  turbulent  flows,  experimentation  gives  0.6  <  PIE  <  1.2.  In  the  present  study, 
PrE  was  set  to  unity,  Pr(-  =  1.0.  The  actual  TELS  operation  is  expected  to  induce  much 
higher  levels  of  turbulence  and  “mixing”  than  are  given  by  the  Schlichting  formula  for 
r-E.  To  simulate  the  high  turbulence  level  induced  by  the  TELS  arm,  pE  was  increased  by 
a  factor  of  10  in  the  numerical  investigations,  and  the  constant,  k^,  was  increased  to 
0.285  for  some  calculations.  The  results  of  these  calculations  are  shown  in  Fig.  6. 


A6.0  MODELING  TECHNIQUES  FOR  BOUNDARY  CONDITIONS 

The  computer  simulation  solves  a  set  of  equations  for  u,  v,  T,  and  P,  which  are  then 
specified  on  the  boundaries  of  the  region.  These  boundaries  include  the  upstream, 
downstream,  free-stream,  and  ground  plane  limits,  plus  the  physical  boundaries  of  the 
engine  and  deflector.  Tire  specific  values  assigned  to  these  variables  are  as  follows: 
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Upstream:  u  =  (crosswind  velocity,  ft/sec) 

v  =  0 

P  =  P  (2116.0  psf ) 

T  =  T  (530°R) 

CO 

Downstream:  Either 

3u/3x  =  0 
3v/ 3x  =  0 
P  =  P 
3T/3x  =  0 

Free  stream:  u  =  U 

OO 

v  =  0 
P  =  P 

CO 

T  =  T 

OO 

Ground  Plane:  u  =  U  (slip-velocity,  neglecting  the 

00 

boundary  layer) 

v  =  0 
3P/9y  =  0 
3T/3y  =  0 

u  =  0 

v  =  400  ft/sec 
P  =  2116.0  psf 
3T/3y  =  0 

Exhaust:  u  =  0 

V  =  VJET  t1’900  ft/sec) 

T  =  1 , SOO^R 
P  =  2116.0  psf 


Engine: 

Inlet: 
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Sides:  u  =  0 

V  =  0 

BT/Bx  =  0 
BP/dx  -  0 

Deflector:  along  top  surface 

u  =  0 
v  =  0 

T  =  1,500°R 
BP/By  =  0 

along  lower  surface 
u  =  ud  (ID,  JD) 
v  =  vd  (ID,  JD) 

T  =  1 ,500°R 
BP/dy  =  0 

The  ability  to  model  deflector  geometry  is  provided  through  specifying  the  velocity 
u<i,  vd,  on  the  lower  surface  of  the  deflector.  Thus,  a  velocity  distribution  is  input  to 
model  the  effects  of  a  physical  deflector.  Methods  exist  for  defining  ud,  vd  profiles  (see, 
for  example.  Ref.  I  7).  In  the  present  study,  the  following  functions  were  defined  for  ud 
and  vd  : 


Case  1 


Case  2 


ud 

(ID,  JD) 

=  0.0 

vd 

(ID,  JD) 

-  0.0 

ud 

(ID,  JD) 

-  -u* 

JET 

sin  ( 

t  7T 

'  2 

ud 

(ID,  JD) 

-  u* 

JET 

sin  ( 

t  IT 

'  2 

vd 

(ID,  JD) 

=  _v* 

JET 

sin  ( 

IT 

v  2 

vd 

(ID,  JD) 

=  V* 

JET 

sin  l 

^  2 

Z~ 

where  the  lengths  2+,  2"  are  local  normalized  lengths  along  the  deflector  surface  (Fig.  6). 
The  velocities  and  V*J1;:T  are  defined  such  that 

V  JET  =  J  U*fET  +  V*JET 
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A7.0  REMARKS  ON  THE  OVERALL  THEORY 

This  numerical  solution  technique  is  a  new  approach  to  solving  the  Navi er- Stokes 
equations  directly  by  the  application  of  existing  theory.  It  has  not  been  fully  tested:  only 
some  preliminary  calculations  have  been  made  for  developing  laminar  boundary  layers 
with  and  without  heat  transfer.  Therefore,  more  research  is  required  to  assess  the 
accuracy  of  this  approach  and  to  incorporate  the  pressure  equation  into  the  evolving 
solutions  to  examine  the  effects  of  compressibility  on  the  jet  deflection  problem.  Some 
standard  flows  -  such  as  parallel  jet  mixing,  turbulent  pipe  flows,  and  diffuscr-type 
flows  -  should  be  computed  and  compared  to  existing  numerical  solutions  and  data. 
Because  of  the  method’s  simplicity  and  broad  applicability,  it  merits  further  investigation. 
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NOMENCLATURE 

ii i  tj,.  U2<j)  Coefficients  m  Ihe  standard  form  of  the  governing  equation  [  Eq.  (  A 2.6 ) ] 

bi  <|,,  l>2 q>  Coefficients  in  the  standard  form  of  the  governing  equation  |  Eq.  (A2.6)) 

Cp  Specific  heat  at  constant  pressure 

DX.DY  Mesh  spacings  in  the  finite-difference  simulation 

dq,  Source  terms  in  governing  equations  (Fq.  f  A2.6)  I 

g  Gravitational  acceleration 

i.  I  Indices  for  mesh  nodes  m  axial  direction 

j,  J  Indices  for  mesh  nodes  m  transverse  direction 

k|  Effective  thermal  conductivity 

k n  Constant  in  viscosity  function 

2+,  2"  Distances  along  deflector  plate  (Fig.  8) 

P  Effective  Prandtl  number 

p  Pressure 

R  Radial  Distance 

Rn  j  Turbojet  exhaust  jet  radius 

T  Temperature 

Frec-stream  velocity 

Li*j  |  T  Maxiinuin  x-direetion  velocity  on  deflector  plate 
u  x-direetion  velocity 

Uj  \-direction  velocity  along  deflector  plate 

V i  [.■  j  Turbojet  exhaust  jet  velocity 
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V* |  ].  t  Maximum  y-direclion  velocity  along  deflector  plate 

v  y-direction  velocity 

vD  y-direetion  velocity  along  deflector  plate 

X  Axial  direction 

y  Transverse  direction 

p  Volume  coefficient  of  thermal  expansion 

/jt|.  Effective  dynamic  viscosity 

i>y  Effective  kinematic  viscosity 

7T  Pi  (3.1415...) 

p  Density 

<t>T  Energy  dissipation  function 
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